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Abstract. - Compressed copper powder has a very large electrical resistance (IMS!), due to 
the oxide layer on grains (100/im). We observe that its voltage-current U-I characteristics are 
nonlinear, and undergo an instability, from an insulating to a conductive state at relatively small 
applied voltages. Current through the powder is then noisy, and the noise has interesting self- 
similar properties, including intermittency and scale invariance. We show that heat dissipation 
plays an essential role in the physics of the system. One piece of evidence is that the instability 
threshold always corresponds to the same Joule dissipated power whatever the applied stress. 
In addition, we observe long-time correlations which suggest that thermal expansion locally 
creates or destroys contacts, and is the driving mechanism behind the instability and noise 
observed in this granular system. 



Introduction. - For over a century pQ, electrical transport in metallic powders has gen- 
erated interest p]. These powders have fascinating properties, such as extreme sensitivity 
to electromagnetic waves, highly nonlinear U-I characteristics, hysteresis, and 1// noise, for 
which fully satisfactory explanations are still lacking. The experiments presented here were 
motivated by the work of Kamarinos et al. p] on compressed copper powders. These authors 
observed an insulating-to-conducting transition at rather low pressure-dependent voltages, 
associated with strong 1// noise. Just above the conduction transition threshold, we ob- 
serve slow temporal evolution of the powder resistance with a noisy component. Electrical 
breakdown of the oxide layers on grains has been invoked for the transition, but this 
explanation is unsatisfactory, for we observe that the noise involves both increasing and de- 
creasing the electrical resistance of the powder. In this letter, we show that this electrical 
noise has interesting scale invariant properties. Scale invariance [5] occurs for various physical 
signals: turbulent velocity financial stock market data [Jj, earthquake energy release 
or world wide information traffic j^j. One of the goals of this study is to use the analytic 
tools developed for turbulent signals in the case of this granular system. We shall focus the 
comparison on two aspects: The statistics of current increments, I(t + r) —I(t), depending on 
the tested time scale r, and the correlations between the amplitudes of these increments. This 
latter will show that the electrical noise in this granular system has a hierarchical organization 
through the time scales. 

(*) E-mail: Eric.Falcon@ens-lyon.fr, URL: |http: / /perso. ens- lyon.fr/eric. f alcon/ 
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Experimental setup. - Two kinds of experiments are performed with commercial copper 
powder samples of 100 /im "spherical particles" P]: U-I characteristics on one hand, and noise 
and relaxation measurements on the other. In both cases, the powder samples are confined 
in a polymethylmethacrylate (PMMA) cylinder, of 10 mm inner diameter, capped with two 
metallic electrodes (stainless steel or brass cylinders). The container is filled with powder up 
to a height of 5 mm, roughly corresponding to 500 000 particles. For the U-I characteristics, 
a sensor measures the force applied to the powder through the electrodes. We occasionally 
embedded two wires inside the powder to check that the resistance is not controlled by the 
electrode-powder interface. Generally, before each experimental run, the container is refilled 
with a new sample of powder. This procedure ensures better reproducibility than simply 
relaxing the confining pressure and shaking the container. 

U-I characteristics. - The DC current, /, is provided by a Kepco Power Supply (BOP 
50-4M). We first apply a static force, F, to a new sample. Then we measure the voltage, U, 
across the sample, as a function of increasing values of the current, /. A single run typically 
lasts 10 s. 




U (V) ~ (Ul) 1/2 (W 1/2 ) 



Fig. 1 Fig. 2 

Fig. 1 - Normalized U-I characteristics for various applied forces: F = 640 (□), 700 (o), 750 (*), 800 

(V) and 850 (A) N. Slope of unity (-). Slope 2 ( ), i.e. I„i/I = 1/2, is close to the instability 

threshold (see text for details). 

Fig. 2 - The relative nonlinear part of the current, I„i/I, versus (UI) 1 ^ 2 . The -^--mark is the point 
where = 1/2, (UI = 0.12 mW ; UI„i = U 2 /R = 0.06 mW), close to the instability threshold 

(see text for details). Symbols are as in fig.^ 

Figure displays the U-I characteristics for various applied forces ranging from F = 640 
to 850 N. / is normalized by the sample resistance Rq(F) at low U, and thus the slope of 
each characteristic is 1 at the origin. We see that at a higher applied force, the departure 
from linear behavior occurs at lower voltages (see fig-GJ. Note that Ro(F) decreases with F, 
1MO being a typical value. We define I n i = I — (U/Rq) as the departure from linearity, that 
is the nonlinear part of /. Figure then shows that plotting I n i/I as a function of (UI) 1 / 2 
collapses all data. At first glance, a simple interpretation for this (UI) 1 / 2 dependence can be 
developed. Indeed, developing I in powers of U such as: 



/ = (U/R ) + cU 2 , yields I nl /I ~ R a cU ~ cR 3 J 2 (UI) 1/2 . 
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However, due to the symmetry of the system, c should be zero. Thus, this interpretation does 
not hold. We prefer to focus on the main observation: The relative nonlinear component of 
the current, I n i/I, seems to depend only on the total dissipated power in the sample, UI. An 
interpretation along these lines will be proposed in what follows. 

Transition and relaxation. - In a new series of experiments, the PMMA cylinder is filled 
with 5 mm of powder, then vibrated and embedded in a cylindrical brass press with 5 mm 
thick walls. The pressure on the sample is generated by means of the lid press acting as a 
screw. At constant low voltage (U ~ 0.5 V), the sample resistance, Ro, is monitored during 
the stepwise pressing of the sample until a maximum pressure, P m , is reached. Rq is found 
to decrease with increasing P m , and to reach a value ranging from 100 kfi to 500 kfi, at the 
end of the pressing process (half an hour later). We then let the system relax for one day. 
After checking that Rq and P m are constant in time, a fixed voltage, U, is suddenly applied 
to the sample, and the current, /, is monitored. If U is smaller than a threshold value U c , the 
sample is in a weakly conducting state with very slow temporal evolutions of /, if any. When 
U > U c , an instability occurs. / rapidly increases at constant U which can be interpreted 
as the resistance relaxing down. The larger U is, the faster the resistance relaxes. Similar 
features have been observed by Kamarinos et al. pj]. 

When the above experiment is repeated for different values of P m , it shows that both U c and 
Ro depend on P m . However, the critical ratio U 2 /R is found to be independent of the applied 
pressure, with a value close to 0.07 mW. This value is close to the point where I n i/I ~ 1/2 
(see -^--mark in fig. El- All these observations suggest that this spontaneous transition from an 
insulating to conducting state is a thermal instability. To further characterize the phenomenon 
observed for U > U c , let us make three remarks: 

• Due to the strong sample pressure, vibrations have no effect; even a strong jolt applied 
to the press has no visible consequence on the signal. 

• We occasionally followed the relaxation down to a resistance of 500 fl, which is almost 
three orders of magnitude lower than the initial one R ; however, this resistance remains 
much higher than a metallic contact would produce between the grains, no matter how 
small the contact [2]. 

• The temporal evolution of the current, at constant U > U c , is not monotonic. Both 
increasing and decreasing events occur for the current. However, the former dominate 
and control the global evolution. This is contradictory to what would result from a 
cascading electrical breakdown for the oxide layers, as supposed by Kamarinos et al. 

Finally, we observe that the direction of the global evolution of the resistance depends 
on the applied voltage. As shown in fig. we first apply to the sample a voltage U = 4 V 
(U 2 /Rq = 0.1 mW > U 2 /Ro), which triggers the relaxation. Half an hour later, we decrease U 
to 0.5 V, and the conductance goes down. Two hours later, we increase U to 3.5 V and we see 
that the conductance goes up again. Consequently, with a well choosen applied voltage, we 
can obtain an approximately constant conductance. Note that these observations are coherent 
with the last item above. We exploit them in what follows. 

Noise. - Our goal is now to obtain a quasi-stationary current signal in order to apply the 
usual tools of signal processing (e.g. spectral analysis), and more sophisticated ones developed 
in the framework of studies of turbulence time series [§]. To start, we consider an initial sample 
resistance of Ro = 0.5 Mf2. To obtain the relaxation, we apply U ~ 7 V (U 2 / Ro — 0.1 mW 
> U 2 /Ro). Five minutes later, the resistance reaches 830 SI, we then apply U = 0.5 V, and 
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- Temporal evolution of the sample conductance. The 4V first applied are above the threshold, 
second part, 0.5V are applied, and finally 3.5V. 



Fig. 4 - Quasi-stationary current noise can be obtained with both increasing and decreasing events. 



the resistance goes up, reaching 1 kO one hour later. Applying now U = 2 V, the resistance 
slowly decreases to a value of 950 after one hour. Finally, applying U = 1.6 V leads to an 
approximately constant current for hours, with stochastic fluctuations (see fig. [3J . At fixed 
voltage U = 1.6 V, the current I is recorded through an acquisition system, with a sampling 
frequency of fn = 25.6 kHz (resp. ft = 128 Hz), the signal being previously filtered at 10 
kHz (resp. 50 Hz) to avoid aliasing. The signal is recorded during 20 s (resp. 65 min) leading 
to a file of 0.5 Mpt. This type of data acquisition was repeated 20 times, first at fjj, then at 
ft, to extract averaged quantities, due to the quasi-stationarity feature of the signal. 

Figure |U displays the log-log power spectra of filtered signals of current recorded at ft and 
Jh- At first sight, this power spectrum of current fluctuations seems to be a power law of the 
frequency. However, when one examines the spectra carefully, a small curvature appears in fig. 
El Letting r designate a time lag or a time scale, we define the i-th order structure function 
Si(r) = ^[/(t + r) — I(t)] 1 ^, where (•) represents an average over time t. We focus on the 
structure function of order four as a function of /j?r, as shown in fig. El In this log-log plot, 
high sampling frequency data present a power law dependence with /jjt. Deviations from 
this are observed at large times (e.g., log 2 /jjt > 16), that is at times r greater than a critical 
time scale t c ~ 3 s. This critical time r c can be understood as a typical effective diffusive 
time of a thermal pertubation within our typical size sample, L ~ 2.5 — 5 mm, estimates 
corresponding to t c ~ 1 — 10 s. This order of magnitude agreement supports the hypothesis 
of a thermally driven phenomenon. With similar reasoning the heat diffusion in a single grain 
provides the short time limit, n n f ~ 0.1 ms, which is on the order of the inverse of our 10 kHz 
filtered frequency. For high frequency data, S±{t) cx r" 4 , with a± ~ 0.57. The same data give 
52 (t) oc t™ 2 , with Q!2 — 0.31. These power-law behaviors are consistent with scale invariance 
over more than 3 decades in time. Intermittency, corresponding to 014 < 2a2, appears more 
clearly from the direct examination of probability density functions, and correlations, in the 
next paragraphs. Here, and 2«2 differ only by one standard error. 

Probability density functions. - We concentrate now on the high frequency part of the 
signal, which presents nice scale invariance properties for r in f < r < t c . Figure shows, 
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Fig. 5 Fig. 6 



Fig. 5 - Average power spectra of current fluctuations sampled at low (V) and high (A) frequencies. 

Fig. 6 - The fourth moment of current differences, Sa(t), versus the nondimensional time scale, Jht. 
Same symbols as fig. El 



for three different r, the probability density functions (PDF) of current differences, 5I T (t) = 
I(t + t) — I(t), normalized to their respective root mean square, a. The exact shape of these 
PDF is rather sensitive to the statistics (see PDF tails in fig. EJ, and would have led to 
smoother curves with a greater quantities of data. However, two remarks can be made in the 
light of what is known for velocity signals in turbulence jH] : 

• The PDF shape changes with the time scale r - this is a direct signature of intermittency, 
as 5*4 cannot be proportional to S% and thus ^ 2a 2 - 

• Current difference PDF are symmetric - time reversal, which a priori should not be 
invoked here, is the only symmetry able to lead to this behavior. Turbulent velocity 
differences have a skewed PDF, S3 being proportional to the dissipated power How- 
ever, even in turbulence, global quantities, like the total dissipated power (equivalent to 
what we measure here) have rather symmetric time difference PDF |10| . 

Multiplicative cascade. - A signal with scale invariance is called self-similar since its 
statistical properties can be described by the same laws at various scales. If the shape of 
time difference PDF changes across scales with a self-similar law of deformation, this signal is 
generally described either in terms of a "multifractal set of singularities", or a "multiplicative 
cascade" across scales, both approaches being considered equivalent [TT] . In the "multiplicative 
cascade", a given scale (e.g., average gradients on a time interval) conditions smaller scales 
(e.g., smaller intervals within this one) in a random Markovian way. It induces not only inter- 
mittency (evolution of distribution shape across scales), but long-range correlations between 
short intervals JlJEl- Therefore, if our current signal can be described by a self-similar 
multiplicative cascade, then the mean squared deviation of ln|W T | should linearly depend on 
lnr, such that 

C 2 (t) = ([hx\8I T \- < \n\5I T \ >f) = -ii\n(r/T) + c, 
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Fig. 7 Fig. 8 



Fig. 7 - Probability density function of the current differences for various time scales: r =0.15, 2.5 
and 40 ms (from top to bottom). Factor of 2 shifts have been applied for clarity. 

Fig. 8 - Comparison between the variance of the logarithms of current differences (C2) and their 
correlation (Lg with fnO = 2 2 ). Both are close and depend logarithmically on the time separation r, 
in agreement with multiplicative cascade models (see text). Solid line of slope — /iln2 = — 1.410 -3 is 
shown. 

where T is some large time scale. Also, as shown by 52, the log-correlation between two 
short intervals (of size 9) should depend linearly on In r, 

L 2 6 (t) = ([\n\6I e (t + T)\\n\5I e (t)\- <\n\5I e \ > 2 ]) 
= -//ln(r/T)+c', 

with the same coefficient /x, c and d being two constants. 

Figure El shows the two experimental quantities C%(t) and Lg(r), shifted by appropriate 
constants. The agreement between them and the linearity in lnr strongly supports a multi- 
plicative cascade description. As in the case of turbulence, large scale events condition those 
at a smaller scale. Any physical interpretation of the phenomenon discussed in this paper 
must address these points. 

Interpretation. - We can now take stock of all of our results and propose a physical 
picture of what is taking place. As shown above, the driving parameter is the total dissipated 
power. This suggests either local heating, able to change the electrical properties of contacts 
(several hundreds of degrees are needed in this case), or thermal expansion, which locally 
creates or destroys contacts. The power involved is of the order of 10~ 4 W. When divided 
by the number of contacts, it remains so small that only electrons can undergo significant 
heating. Such "hot" electron effects have been reported in systems having some analogies 
with this one However, the influence of the large scales on the small ones, as shown by 
the observed logarithmic correlations, cannot be taken into account by such a local process. 
Therefore, thermal expansion seems to be the mechanism driving the instability and the 
associated noise. 

In this spirit, we attempt to show the correspondence between the formal multiplicative 
cascade and electrical conduction in our powder. Since the contact distribution in a powder 
is very inhomogeneous, one would also expect an inhomogeneous current distribution. Thus 
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Joule heating should create inhomogeneous increasing stresses in the powder. Very small 
thermal expansion can result in dramatic changes in the current paths, thus in the distribution 
of this Joule heating, and so on. Such events can occur at any scale, ranging from the size 
of the sample and the grain size. The large scale events should influence the small ones, as 
suggested by our study. A full confirmation of this idea requires longer carefully controlled 
studies. 

Conclusion. - The results of this work are two-fold. First, we show that the sponta- 
neous decrease in resistance of a copper powder sample above a voltage threshold is due to a 
thermal instability, and not to electrical breakdown, as previously had been proposed. This 
conclusion results from the observation that dissipated power drives the phenomenon, in spite 
of the probable smallness of the induced temperature inhomogeneities. Second, we propose 
a procedure yielding an interesting self-similar process in this non-equilibrium system. Our 
system displays both intermittency and multiplicative cascade-like two point correlations, in 
ways that are interesting to compare and contrast with the archetypical case of turbulence. 

* * * 
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